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ABSTRACT 

In a Brans-Dicke (BD) cosmological model, the energy density associated with 

some scalar field decreases as a V / with the scale factor a(t) of the Universe, 
giving a matter with an Equation of state p = — f — — — - ] p. In this model, the 

3 \l+L} J 

Universe could be closed but still have a nonrelativistic-matter density corresponding 
to its critical value, o = 1. Different cosmological expressions, such as, luminosity 
distance, angular diameter, number count and ratio of the redshift tickness-angular 
size, are determined in terms of the redshift for this model. 
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1 INTRODUCTION 

In the Friedman-Robertson- Walker (FRW) model, there are 
enough evidences that the Universe is flat, with a large com- 



ponent lot negative pressure. This component was consider 
at first to be just Hits coamological consianl (or vacuum en- 

ergy) (Tuner et al. 1984; Peebles 1984; Ostriker & Steinhard 
1995; Liddle et al. 1996). Another possibility was to consider 
topological deflects (Vilenkin 1984), but, the one which has 
receiving a great deal of attention today, is related to a scalar 
field, Q, the so-called "quintessence" model or the QCDM 
model (Caldwell, Dave & Steinhardt 1998). This scalar field 
is characterized by a very negative pressure, i.e., Pq = wpQ, 
with w < —1/3. 

Measurements of distant SNe la (at z ~ 1) indicate that 
the expansion of the Universe is in accelaration rather than 
deceleration (Perlmutter et al. 1998). 

According to Garnavich et al. (1998) this is consistent 
with the existence of an unknown component for the en- 
ergy density, which could be considered to correspond to 
the quintessence component. On the other hand, a test of 
the standard model, including spacetime geometry, galaxy 
peculiar velocities, structure formation, and early Universe 
physics, supports, in many of these flat Universe 

model with the presence of a cosmological constant (Pee- 
bles 1998). Specially, a model in which a short period of 
inflation occurs at a very early epoch in the evolution of the 
Universe. Most of these models predict that the total density 
parameter, Sl , be unity. 



Given the idea that the Universe could be described by 
a flat geometry, an interesting question to ask is whether 
this flatness could be due to a local effect. T his sort of ques- 
tion has been considered in the literature ((Kamionkowski 



& Toumbas, 1996)). There, a k = 1 model was taken into ac- 



count, together with a total density parameter correspond- 
ing to il < 1. In this model, openness is obtained by 
considering a matter component whose equation of state is 

p = — —p. This sort of state of Equation has been reported in 

models in which topological defects, such as texture or tan- 
gled strings, are important components of the total energy 
of the Universe. 

In this paper we analyse a FRW closed model, k = 1, 
filled with dust and a effective density energy component 
characterized b y a negative pressure, using the Brans-Dicke 
(BD) theory ((Brans & Dicke, 1961)) with a potential as- 



sociated to the BD field. In particular, we investigate the 
conditions in order to have a model with positive curvature 
which mimics a flat universe at low redshift. This imply 
to determine the contributions of the scalar field Q(t) and 
the BD field to the total energy density, which cancel the 
contribution due to curvature fl k = —1/aoHo. We should 
note that the obtained model is the generalization of the 
Einstein-de Sitter model (k = 0) in the context of the BD 
theory. Our model is far from being realistic, since it gives 
an age for the universe,io, very close to 2Ho/3, and a decel- 
eration parameter, qo, very close to 1/2. These values, as we 
mentioned above, desagree with the measurements of dis- 
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tant supernovae. Nevertheless, we clearly pointed out that 
our proposal is to set the conditions under a closed uni- 
verse mimics a flat universe in the BD theory. At the end 
of the following section we shall briefly discuss how in the 
BD theory a flat universe with negative pressures, presents 
acceleration in agreement with actual observations. 

We studied this sort of model in a previous paper 



((Cruz, del Campo & Herrera, 1998)), but that analysis was 



limited to the model's intrinsic characteristics, such that the 
explicit determination of the scalar fields, Q(t), its potential 
V(Q) and the BD potential V($), where $ represents the 
BD scalar field. In this paper we investigate the model's cos- 
mological characteristics, such that, the proper distance to 
the horizon, the luminosity distance, the angular size, the 
differential number of galaxies, and the ratio (l/z)5z/58, 
where, all of these quantities are given as a function of the 
redshift z. We compare these parameters with that corre- 
sponding to the Einstein-de Sitter model. 

We should note that our model, in the limit uiq — » oo 
and $ = const, gives rise to the k — Einstein-de Sitter 
model. 



2 CHARACTERISTICS OF THE MODEL 

Assuming homogeneity and isotropy, the FRW metric for a 
closed universe is 

ds 2 = dt 2 — a(t) 2 dSii—i, (1) 

with df2j. =1 representing the spatial line element asociated to 
the hypersurfaces of homogeneity, corresponding to a three 
sphere. a(t) represents the scale factor, which together with 
the assumption that the Q scalar field is homogeneous, i.e., 
Q = Q(t), we obtain the fundamental field Equations of the 
BD model, given by (where the dots representing derivatives 
with respect to time t. We use units in which c = % = 1 ) 




8t , 1 U(&) 



and 

$ +3H $ 



Stt 



jPM, 



(6) 



(7) 



2lo + 3 

respectivelly. Equation (jij) remains unaltered. Note that this 
set of equations mimics a flat (fi m = f?o = 1) Universe. 

Assuming that the matter content, pM, is dominated by 
dust, i.e. pm ~ a -3 , the solutions of Equations ([]) and (Q) 
are the well known power law solutions of a flat BD model, 
given by 



a=a °(lY 



and 



$(t) = $ 



(8) 



(9) 



where t is the current age of the Universe. Using the usual 
formula for the redshift z, 1 + z = a a /a, we find that 



t 



2 ( LJ + 1 

3 I Wo + 4/3 



H^il + z) 



(10) 



where H = yj Sixp M /3<^ , is the current Hubble constant 
taken as H = WOhkms -1 Mpc~ x , with h in the range 
0.6 < h < 0.8 (Bureau, Mould & Staveley-Smith 1996). 
From Equation (llOj) we can obtain the age of the Universe, 

to = I (^TITs) H ° X = (zJtTITs) tE ° ' where if is the 

Einstein-de Sitter value for the age of the Universe, given by 
. At first glance, the factor appearing in t seems to 



LOo + 1 



< 1. 



decrease the age of the Universe, due to . 

\lj + 4/3 / 

But, since u> > 500 (Reasenberg et al. 1979), this factor will 
be small, and almost equal to one, and, therefore, we could 
use to — to . 

The growth of the Q-field is given by 



Q J - 2Wo Uo + l 



(11) 



$ +3H $ 



and 



$ 3 d 



2lu + 3 d$ ^ <E> 2 J 

8tt 
2w D +3 



+3H Q= 



dV(Q) 



[pM + (1 - 3w) p Q ] , (3) 



(4) 



The condition describing a model mimicing a flat Uni- 
verse is (Cruz, del Campo & Herrera 1998) given by 



d (U{$) 



1 [/($) 



d$ ^ $ 2 j 2 (Uo + 1) $ 3 

+ i!(TT^)(^) 

where $o is the the actual value of the BD scalar field. Under 
these conditions Equation (H) and Equation (H)reduce to 



and its present value is then determined by the values of the 
Hubble constant, H a , and the BD-parameter uj . At large 

values of u) (lu — > oo) it becomes (Q /0) = 2H . On 

the other hand, the potential V(Q) associated with this field 
is given by 



V(Q) = M 



QJ 



where M is a free parameter. 

„. 3 / uj + 1 

Since po = - — ■ 

^ 2 U„ + 5/4 



parameter S7q° is defined by Q,q = 



(12) 

V(Q) and the actual density 

Z*PQ 



3$ m ' 

o; + 1 \ 47rM 4 / w + 1 



we find that 



SIq, where Qq 



k (j + 5/4 J QoH'i \uj + 5/4 / 
corresponds to the scalar field density parameter in Einstein- 
de Sitter's model of the universe. 
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If we consider the same range for J7q in our model, we 
find that this imposes a constrain on M. Since today it is 
required that V(Q ~ M p i) » p m ~ 10~ 47 GeV 4 , where M p * 
is the Planck mass, we obtain 



M 



iM, 



(13) 



This gives M ~ 10 Gel/ for oj ~ 500. This value is compa- 
rable to particle physics scales. We mention that the evolu- 
tion of Q and the values of Qq today are very insensitive to 
the initial value of Q, due to its attractor solution (Zlater, 
Wang & Steinhardt 1998). 

One point that needs to be taken into account in this 
kind of model is the fraction of the Universe in causal con- 
tact. To do so, we employ a comoving observer at coordi- 
nate (r = 0, 8, cp) at time t. A light signal satisfies the 
geodesic equation of motion da 2 = 0. Therefore, a light sig- 
nal emitted from (ri,6,<p) at time t = 0, following the line 
6 = 6 — const., will reach the observer at time t 



dt' 



dr 



O(t') J (l_ r 2)3 
o o 

and because the proper distance to the horizon is 



(14) 




Redshift z 



d H (t) = J (g rr ) 2 dr, 

o 

it is found that 

dH(z)=a ° (1+,) 

where a (w , H 



-|8(« )1 



(15) 



(16) 



do H 



+ 3/2) (a; + 4/3)/ (w + 2) 

and /3 (d ) = — ( UJ °~ >r '^ \ j n pjg ^ we nave plotted du{z) 
2 V 0-? o + 1 / 

as an function of z. In Einstein-de Sitter model we have 
nearly obtained the same curve. 

Since, uj 3> 1, we can compare du{z) in the BD the- 
ory with the corresponding expression in Einstein's theory 
of general relativity (Einstein-de Sitter model). Expanding 
Equation (jl6]) up to the first order in l/ui , yields 

: d%(z) + g{z,H )(l/u„) + O{l/u )\ 



du{z) 

where g(z, H ) becomes defined by 



g(z,H ) 



h: 



(l + z) 



VTT. 



1 + 



£n(l + z) 



(17) 



(18) 



and dff(z) represents the horizon distance in the Einstein- 
de Sitter model. This expression presents a maximum for 
2 ~ 2, and there, the difference, Adtf (z) = dij(z) — d^{z), 
computed up to order 1/lj , becomes a maximum. Its value 
is not significant, since it is less than one percent. 

In the following we determine the deceleration param- 
eter q for our model. This parameter is defined by q = 
— (a/a// 2 ) . By using the solution given by Eq. (|^), we ob- 
tain that 



_ 1 fiOo + 2 



(19) 



Figure 1. This plot shows the horizon distance dn(z) (in unit 
of H ) as a function of the redshift z, for the BD parameter 
lu = 500. The curve coincides with the corresponding to the 
Einstein-de Sitter model. 



Note that we can write q — 7 (w ) — 1, where 7 (w„) is the 
inverse of the exponent in the expression of the scale factor 
3 / u + 4/3 s 
2 ^ lu + 1 

bound for cu , i.e. cu — 500, we find that the deceleration 



a(t), i.e. 7 (a> ) 



If we consider the lower 



parameter has a value close to one half, as it should be in 
the Einstein-de Sitter model 

Following the approach done by Uehara and Kim 
(1982), we may write directly an expression for the decel- 
eration parameter q given by 



<7o 



. u 2 
e + Y e 



uJo + 3 
loo + 3/2 



(20) 



Notice that this factor is related to the present rate of 
change of the Newton's gravitational constant expressed by 

(G/G) ( =| J |o, since G($) = ~. 



$ !U ' v ' $ 

The contribution to the deceleration parameter in 

Equation (^5|) is small, since its experimental upper limit 

is given by /$J < 10~ 1( V _1 (Helling et al. 1983; 

Dickey, Newhall & Williams 1989; Shapiro 1989). More re- 
cent, measurements on white dwarfs, have decreased the up- 
G 



per limit to 



G 



!o< 10" 



yr 



(Garcfa-Berro et al. 1995). 



Still smaller upper limits for this quantity have been re- 
ported (Miiler et al. 1991), where lunar laser-ranging stud- 

G 

ies of the moon's Earth orbit yields | — |o< 10" yr~ . 
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Since eo 1, we can write, q 



uJo + 3 



fi TO , which 



2 + 3/2 / 

becomes exactly one half the result of Einstein-de Sitter's 
theory. 

At this point we would like to mention that it is possible 
to describe in the same spirit, the situation in which the 
model represents an accelerating flat universe, i. e. a model 
in Brans-Dicke theory in which now Oq = fl m + Six = 1. 
But, as we will see, this case becomes quite complicate to 
handle. We shall postpone the details of these studies for the 
near future. We shall restrict ourselves here to give a brief 
description of this situation. In this case the condition Q 
becomes 



di> 



$2 



3w) 



[/($) 



3(l- w )_+3(l 



3w) 



$ 2 a 2 ($) 



(21) 



If we want to get an explicit form of the Brans-Dicke scalar 
potential, U($), we need to know the scale factor a as a 
function of the Brans-Dicke scalar field In order to do 
this, we consider the set of basic field equations 



H 2 + H U /$) = f£ U /$) 



8tt 



and 

$ +3H $= 



8tt 



2lu + 3 



PAI + 



2A 



2lo + 3 



:4>. 



(22) 



(23) 



(together with equation (^)) that can be solved exactly (Ue- 
hara and Kim 1982). The solution for A > is given by 



o(t) 



4tt 

T 



A cosh (jjAi) 
Btanh^At) - y / (4 7 r/A) 2 - A 2 



A 2 



Btanh (^Ai) + x /(4tt/A) 
where, ry 2 , A 2 , B, a(u) ) and 0(u>o) are given by 



, (24) 
2(4 + 3lo ) 



3 + 2u; 



4tt 
T 



47T 



2A 4 + 3w 

1 + LOo 



9o 

Po 



and 

4 + 3w 
The interval of time, 



[(l+u )6 + H } 2 {> 0) 
1 / 3 + 2lu 



4 + 3lj d 
At = t 



tc 



, respectively, 
is related to 



tc = 



4tt/A - A 



f±_A 



Expressions ( |24[ ) together with equation ( |22[ ) allow us 
to determine the scalar field $ as a function of time. Thus, 
in principle, we could get the scale factor as a function of 
the Brans-Dicke scalar field. 

We should note that the deceleration rate q becomes 
in this case 



<lo 



Up 2 , 1 I Up + 3 

£o+ 3 e ° + 2 l^ +3/2 
A 



2u) 
2wp + 3~ 



-Six 



(25) 



with Q\ defined by „ . Notice that if the cosmological 

O-Oo 

constant contribution to the total density parameter is sig- 
nificant, in agreement with the Supernova observations, the 
deceleration parameter becomes negative and thus the uni- 
verse will show an acceleration. We stop here this analysis 



and, as we have mentioned above, we shall describe in more 
details this situation in near future. 

In going back to our simple Einstein-de Sitter-like model 
in BD theory, we discuss in the following some kinematical 
properties. 



3 KINEMATICS OF THE MODEL 

3.1 Luminosity Distance-Redshift 

The "luminosity distance" is defined as the ratio of the emit- 
ted energy per unit time, C, and the energy received per unit 
time T 

dl = 4^ ( 26 > 
In the absence of an expansion, the luminosity distance is 
simply the physical distance to the source. In closed FRM 
cosmology, the luminosity distance to a source at coordinate 
(ri,9, <j>) is, assuming for convenience that the observer is 
located at r = 0, 

d\(z) =alrl{z){l + zf, (27) 

The factor a^r 2 in this expression is nothing but the "inverse 
square law" , since a two-sphere surrounding the source en- 
compassing the observer has an area equal to 4-Ka 2 r 2 (z). The 
factor (1 + z) 2 appears due to the redshift of the radiation 
between the time of emission and the time of detection. The 
parameter n is determined by the expression 



ri(z) 



dr 



l 



dz 



where E{z) is given by 



E{z 



(28) 



(29) 



a H j3 {u) ) a (uo , H a ) ' 

Thus, we find for the Luminosity distance-redshift rela- 
tion 

d L (z) = a (z + 1) sin {a (w , H a ) [l - (1 + z)- 0M ] } . (30) 

Fig. 2 shows how the luminosity distance, dh (z), changes 
with redshift parameter z, for Einstein-de Sitter model (dot- 
ted line) and BD closed model (solid line) . From this Fig- 
ure we observe that the BD-theory begins to differ from 
Einstein-de Sitter's theory at z ~ 100. If we take the value 
for the redshift z corresponding to the value associated to 
the last scattering surface, i.e. z = zls — 1.100, the lumi- 
nosity distance in Einstein-de Sitter's theory differs approx- 
imately from its analogous expression in the BD theory by 
a few percentage points. 

Note that for small z, Equation ( pic] ) yields 

H d L (z) = z +l(l-p (u„)) z 2 + o(z 3 ), (31) 

where we have rescaled the luminosity distance by a factor 
of a Hof3 (w ) a (w , Ho) . The first term of this expression 
is nothing but the Hubble law, and from the second term, 

we can read the q parameter, q = (3{uJ ) = - ( — ), 

2 \ uj + 1 / 

which coincides with the result obtained above. 

Also, from Equation ( j3o| ) we can use the apparent bolo- 
metric magnitude m(z) of a standard candle (with M equal 
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approaching asymptotically a limiting value. At ui = 500 
this difference yields the value Am ~ 3.32. 

On the other hand, if we consider the values specified 
by Goobar & Perlmutter (1995) in which the apparent mag- 
nitude (R-band) was ma — 22.17 ± 0.05 at redshift z = 0.5 
and rriR = 25.20 ± 0.05 at redshift z — 1.0, we obtain 
Amu = 3.03 ± 0.05. This value represents a similar order of 
magnitude as that obtained for u) = 500. 




) 600 

Redshift z 



Figure 2. Plot of the luminosity distance di J (z) (in unit of H ) 
as a function of the redshift z. The continuous line corresponds 
with u a = 500 (BD theory) and the dotted line represents the 
luminosity distance in the Einstein-de Sitter model. 



to the absolute bolometric magnitude) as a function of the 
redshifts 

m(z) = M + 5log[d L {z)] + 25. (32) 
H cIl(z), then we obtain, from the 



If we define Dl{z 
latter Equation 

m(z) = M +5log[D L (z)}, 



(33) 



where M. = M — 5log(H ) + 25, is the magnitude "zero 
point" , which can be determined observationally. Hamuy et 
al. (1996) have determined the value M. = — 3.17±0.03 using 
18 supernovae discovered by the Calan/Tololo researchers. 
This value is supposed to be independent of the redshift z. 
Therefore, if we take two different values for the redshift, we 
can obtain a variation of the apparent bolometric magni- 
tude which is given by Am(z2,zx) = 51og [dt,(z2)/di,(,Zi)], 
or explicitly 



Am (02, zi ) = 5 log 

J (1 + 22) sin [a (lu , H ) (l- (l + 2 2 ) 



-/3("o) 



)] 



(1 + 21) sin [a {uj ,H ) (l - (1 



+ 21 



.(34) 



This expression could be used to restrict the value of the 
BD parameter u> - 

By using the values 21 = 0.5 and 22 = 1.0, we find the 
Am increases monotonically with the BD-parameter oj . 
When this parameter reaches a value close to u ~ 20, the 
growth of Am with respect to uj increases, but slowly now, 



3.2 The Angular diameter-redshift 

The angular diameter distance dA between a source at red- 
shift 2 2 and 21 < 22 is defined by 



dA (21, 22) = 



a a sin [Ax (21,22) 



(35) 



(I + 22) 

where Ax (21, 22) is the polar-coordinate distance between 
a source at 21 and another at 22, in the same line of sight 

A x (2i, 22) = a (u>o,H ) [(1 + 2i)-" ( "° ) - (1 + 2 2 )"' 3( "° ) ] (36) 

The corresponding angular size of an object of proper 
length £ at a redshift 2 results in 8 (2) ~ t/dA (0, 2), which 
becomes (in units of £H ) 



Q(z) 



(1 + *) 



^oH sin { Q ( Wo> Ho ) [1 _ (1 + Z )]- 0M } 



(37) 



Fig. 3 shows the plot of O as a function of 2 in Einstein-de 
Sitter's theory and the BD theory with u — 500. These two 
theories coincide in the < 2 ~ zls- At z = zls, it is found 
that in the first approximation of 1/uj , their difference is 
AO = Qbd — Qe ~ 147(l/w ), for any large value of 0J o - 



In this plot, we h ave taken aoHo = U 2/5 ( Kamionkowski 
fc Toumbas, 1996[). 



3.3 The number count-redshift 

The number of galaxies in a comoving volume element in an 
angular solid area d£l with redshift between 2 and 2 + dz is 
sensitive to the number of galaxies, n, in a comoving volume 
element dV c and the spatial curvature: 



dNaai = ndVc = n- 



-drdQ 



(38) 



/.From this relation we can write the differential number of 
galaxies per steradian per unit redshift, 



dN ga 



f(*) = 



[Ax (2 



dzdtt y ' H E{z) 

Using the expression for A%(z) from Equation ^ and E{ 
from Equation (|2{]), we obtain 

^T%y (2) = n a 3 a{uj , H )(3(u ) x 



(39) 



{a{w ,H ) [!-(! + 2)-^"°']} 



(40) 



In Fig. 4 we have plotted the number-redshift relation (given 
by Equation ([$9|)), for cu = 500. Here, we have confirmed 
that, for a large range of the redshift z, the function of the 
number of galaxies per unit of steradian per unit redshift in 
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0,01 0,1 1 10 

Redshift z 



1,0 




0,0 -| 1 1 1 1 1 1 1 

0,0 0,5 1,0 1,5 2,0 

Redshift z 



Figure 3. This plot show the angular size (in unit ofJHo) as a 
function of the redshift z, as expressed by Equation fl37| ) for the 
case u> = 500. We should mention that the corresponding curve 
for the Einstein-de Sitter model is indistinguishable in the range 
< z < z LS . 



Figure 4. The differential number of galaxies per unit redshift 
per steradian in units of nH^ 3 for the BD parameter w = 500. 
At a low redshifts these differential numbers rapidly approach 



Einstein-de Sitter's model becomes almost indistinguishable 
from its similar expression in the BD-theory. 
For small redshift one obtaines 



0) ~ n ( Oo a(u , Ho) (3{u )f z 2 + O (z 3 ) 



(41) 



dN gal 
dz dQ. 

This expression indicates that, for small redshift, the slope in 

the plot — — v/s z 2 could provide some information about 
dz ail 

the lu BD parameter, assuming that the present value of the 
Hubble constant H is known. This, in principle, could be 
checked by the corresponding observations. 



3.4 The ratio of the redshift thickness-angular size 

The redshift thickness Sz and the angular size 86 of a spher- 
ical structure that grows with the expansion of the Universe 
have a dependence on z given by 

E{z)sm [A X (z)\ 



1 &z ( \ u 
--(z)=a oHo 



or explicitly, 
1 Sz, 



(42) 



(1 + z) lM sin 



(43) 



_ {a(co ,H )[l^(l + zf { ^]} 

z 89 a H a(uj ,H ) /3(u )z 

1 8z 

Fig. 5 (lower panel) shows how — changes with redshift z 
in the BD-theory, for lu — 500. We see that, in the range of 
< z < 10, there is no difference when it is compared with 



its analogous expression in Einstein-de Sitter model. How- 
ever, at large redshift i.e. z = zls ~ 1.100, their difference 
becomes significant, as we can see from Fig. 5 (upper panel). 
Thus, the two expressions differ at large redshift values. 



4 CONCLUSIONS 

We have studied a model in which the "quintessence" com- 
ponent was included. In this model we have chosen a Uni- 
verse with closed topology (k — 1) in the BD theory, and we 
have added a scalar potential for the BD-scalar field. The 
quintessence contribution to the matter component was fine 
tuned to exactly cancel the curvature term (with the help 
of the BD scalar potential) in the corresponding field Equa- 
tions. 

Under these conditions, different cosmological expres- 
sions were calculated as a function of the redshift z. For 
instance, the luminosity distant, the angular diameter, the 
number count and ratio of the redshift thickness-angular 
size. All these quantities were compared with their analo- 
gous expression related to Einstein-de Sitter model which 
seems flat at low redshift. In some cases (Luminosity dis- 
tance and the ratio of the redshift thickness-angular size), 
the corresponding expressions become distinguishable at a 
high enough redshift. At that redshift, these differences are 
difficult to be directly detected . Perhaps, if we consider 
other observational facts, such as, the anisotropy of the cos- 
mic microwave background radiation, we might say some- 
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Redshifts z 



Figure 5. In the lower panel we plot (8z/S8) /z as a function of 
z in BD theory (for ui = 500) at low redshift. In upper panel we 
have plotted (8z/88) /z in Einstein-de Sitter (dashed line) and BD 
theories (continous line) for z > 10. Note that these two curves 
become different at very high redshift. 



thing about these differences, and use these results for es- 
tablishing a limit for the BD parameter, ui , since the two 
theories differ in the value of this parameter. 
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